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PREFACE 


The  asymptotic  development  of  the  two-dimensional  incompressible 
wake  far  behind  a  flat  plate  is  a  classical  problem  in  fluid  mechan¬ 
ics,  but  it  is  only  recently  that  the  correct  asymptotic  series  solu¬ 
tion  has  been  given.  This  Memorandum  treats  the  analogous  problem 
for  axisymmetric  flow,  which  corresponds  to  the  wake  flow  behind  a 
very  slender  cylinder  or  "needle";  however,  the  solution  has  some  rel¬ 
evance  to  the  wake  behind  any  finite  axisymmetric  body.  The  results 
presented  here  should  be  of  interest  to  those  working  in  the  general 
area  of  wake  flows. 


summary 


This  Memorandum  considers  the  asymptotic  development  of  the  axi- 
symmetric  incompressible  wake  far  behind  a  long,  thin  cylinder.  This 
is  the  axisymmetric  analogue  of  the  classical  two-dimensional  flat- 
plate  problem;  however,  the  resulting  solution  is  applicable  to  a  wider 
class  of  axisymmetric  bodies. 

In  Section  II  an  asymptotic  expansion  in  inverse  powers  of  the 
axial  distance  downstream  of  the  cylinder  is  attempted;  such  an  expan¬ 
sion  is  suggested  by  the  Oseen  solution  to  the  problem.  The  solution 
is  calculated  to  the  second-order  term.  However,  the  condition  that 
the  perturbation  velocity  go  to  zero  exponentially  far  from  the  axis 
is  not  satisfied  by  the  second-order  term.  Hence,  this  term  must  be 
rej  ected  on  physical  grounds. 

A  similar  difficulty  arises  in  the  two-dimensional  flat-plate  prob¬ 
lem  in  the  third-order  term.  Resolution  of  this  difficulty  is  accomp¬ 
lished  by  adding  a  logarithmic  term.  The  following  term  contains  an 
indeterminate  numerical  factor;  this  indeterminacy  can  be  shown  to  be 
associated  with  the  neglect  of  the  initial  velocity  profile  in  obtain¬ 
ing  the  asymptotic  far-wake  solution.  Indeterminate  factors  also  ap¬ 
pear  in  various  other  terms  in  the  series  expansion,  and  these  can  be 
shown  to  be  connected  with  the  existence  of  eigensolutions  to  the  equa¬ 
tions  satisfied  by  individual  terms  of  the  expansion. 

In  Section  III  the  axisymnetric  problem  is  reformulated,  and  it 
is  shown  that  as  in  the  two-dimensional  problem  the  cause  of  the  break¬ 
down  of  the  solution  given  in  Section  II  is  the  existence  of  a  loga¬ 
rithmic  term  in  the  second-order  approximation.  The  logarithmic  term 
is  calculated  explicitly,  but  again,  the  next  term  is  indeterminate  be¬ 
cause  of  neglect  of  the  initial  velocity  profile.  The  next  approxima- 
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tion  is  shown  to  involve  a  (in  x)  term  which  must  be  included  for  the 
next  term  to  have  the  correct  behavior  far  from  the  axis.  Indetermi¬ 
nate  factors  associated  with  the  existence  of  eigensolutions  also  ap¬ 
pear  in  succeeding  terms. 
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I.  INTRODUCTION 


At  the  present  time  the  complete  analytic  solution  for  the  wake 
behind  a  body,  slender  or  blunt,  is  not  obtainable  because  of  the  es¬ 
sential  nonlinearity  of  the  problem  and  the  transition  from  laminar 
to  turbulent  flow  at  some  point  in  the  wake.  In  a  search  for  a  more 
tractable  problem,  one  is  led  to  a  consideration  of  incompressible 
laminar  wakes.  However,  even  for  this  restricted  problem,  an  analytic 
solution  valid  in  the  entire  wake  is  not  possible.  One  then  looks 
further  to  see  if  there  exist  portions  of  such  flows  where  an  analytic 
treatment  is  possible.  There  do,  in  fact,  exist  two  such  regions  for 
certain  bodies:  the  first  of  these  is  the  near  wake  immediately  be¬ 
hind  a  flat  plate  or  very  slender  cylinder,  wuile  the  second  is  the 
wake  far  behind  either  of  these  bodies. 

For  the  flat  plate,  no  recirculating  zone  exists  beyond  the  trail¬ 
ing  edge;  for  the  slender  cylinder,  whatever  small  recirculating  zone 
is  present  may  often  be  neglected  (see  Ref.  1).  It  is  then  possible 

to  obtain  a  series  solution  for  small  distances  from  the  rearmost  point 

(2) 

of  the  body.  The  flat -plate  near  wake  was  first  treated  by  Goldstein, 
while  the  solution  for  the  very  slender  cylinder,  or  "needle,"  was  given 
recently  by  Viviand  and  Berger. ^ 

(3) 

The  wake  far  behind  a  flat  plate  has  been  analyzed  by  Tollmien, 
Goldstein,  ^  and  Stewartson.  ^  This  Memorandum  discusses  the  far 
wake  behind  the  "needle."  (For  reasons  that  will  be  indicated  later, 
however,  the  solution  obtained  will,  in  fact,  be  valid  for  the  far  wake 
behind  any  finite,  closed  axisymmetric  body.)  The  problem  will  be 
treated  by  seeking  an  asymptotic  series  solution.  In  Section  II,  we 
will  attempt  a  series  solution  in  powers  of  x  n,  n  taking  on  successive 

integral  values;  but  we  will  also  show  that  the  second  approximation, 
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the  x  term,  does  not  exhibit  the  correct  behavior  far  from  the  wake 
axis,  and  hence  the  expansion  breaks  down. 

In  Section  III,  following  Stewartson 's  approach  to  the  two-dimen¬ 
sional  case,^  the  problem  is  reformulated  and  the  solution  expressed 
in  a  different  form,  enabling  us  to  explain  the  breakdown  of  the  direct 
expansion  of  Section  II  and,  in  addition,  to  indicate  the  form  of  the 


correct  expansion.  The  second  approximation,  which  was  rejected  in 

_2 

Section  II,  is  given  explicitly;  the  x  term  in  the  expansion  contains 
an  undetermined  constant  connected  in  some  way  with  the  initial  profile 
of  u,  which  is  neglected  in  the  asymptotic  development.  The  nature  of 
the  indeterminacy  in  succeeding  terms,  due  to  the  neglect  of  the  ini¬ 
tial  profile  and  its  relationship  to  the  existence  of  eigenfunction  so¬ 
lutions,  is  also  indicated. 

For  the  two-dimensional  flat-plate  problem,  the  Oseen  approxima- 

(3) 

tion  for  the  velocity  far  downstream  was  first  given  by  Tollmien. 

(4)  _1 

Goldstein  found  the  next  approximation  (of  0(x  )),  then  joined  the 

(2) 

far-wake  solution  to  his  near-wake  solution  by  a  graphical  method. 

In  this  way  he  determined  the  origin  on  the  x-axis  (initially  arbitrary) 

in  the  far-wake  solution.  Goldstein  also  tried  calculating  the  next 

-3/2 

approximation  (of  0(x  ))  but  had  to  reject  his  solution  on  physical 

grounds  because  it  did  not  have  the  correct  exponential  decay  at  infin¬ 
ity.  Stewartson^  showed  that  this  difficulty  could  be  resolved  by 

-3/2 

adding  an  extra  term  0(x  in  x)  into  the  expansion;  he  also  found 

-3/2 

that  a  numerical  factor  occurring  in  the  term  0(x  )  remained  inde¬ 

terminate  and  that  this  indeterminacy  is  probably  connected  with  the 
neglect  of  the  initial  velocity  profile  in  the  asymptotic  solution. 
Stewartson  also  clearly  discussed  the  existence  and  significance  of 
eigensolutions  in  asymptotic  expansions  in  boundary- layer  problems. 

There  seems  very  little  possibility  in  the  axisymmetric  case  of 
reasonably  determining  the  origin  for  x  by  numerically  matching  the 
"needle"  near-wake  solution  of  Viviand  and  Berger^  with  the  asymp¬ 
totic  far-wake  solution,  as  the  near-wake  solution  for  the  "needle" 
does  not  extend  far  enough  downstream  for  such  a  matching  to  be  accu¬ 
rate.  In  particular,  Viviand  and  Berger  find  that  the  range  of  valid¬ 
ity  in  the  downstream  direction  of  the  new  wake  solution  is  of  the  or¬ 
der  of  the  cylinder  radius  and  that  for  certain  ranges  of  the  parame¬ 
ters,  it  is  much  less  than  this  radius.  On  the  other  hand,  for  the 
two-dimensional  near  wake  behind  a  flat  plate,  the  range  of  validity 
of  the  expansion  used  by  Goldstein  is  an  appreciable  fraction  of  the 
length  of  the  plate. 

The  analyses  referred  to  above,  as  well  as  the  present  one,  are 
boundary- layer  analyses;  that  is,  the  usual  Reynolds  number  is  assumed 
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to  be  large  enough  so  that  boundary- layer  approximations  are  valid  in 
the  far-wake  region.  The  more  general  problem,  finding  solutions  of 
the  full  Navier-Stokes  equations  at  large  distances  from  a  finite  body, 
has  been  treated  by  Chang^^  for  two-dimensional  flow  and  by  Childress^  ^ 
for  axially  symmetric  flows  and  general  three-dimensional  flows.  In 
these  analyses,  which  are  valid  for  a  fixed  value  of  the  Reynolds  num¬ 
ber  which  is  assumed  neither  large  nor  small,  a  small  extraneous  param¬ 
eter  is  introduced;  the  construction  of  the  asymptotic  expansion  is 
then  recast  as  a  perturbation  for  small  values  of  this  parameter.  Due 
to  the  presence  of  the  viscous  wake,  the  perturbation  is  in  general  a 

singular  one  and  is  treated  by  the  method  of  "inner  and  outer  expansions" 

(81 

(or  "matched  asymptotic  expansions"'  ')•  The  results  of  the  present 
analysis  agree  with  those  of  Childress  for  the  axisyrametric  case,  but 
there  are  significant  differences  in  the  basic  approach,  details  of 
the  analysis,  and  interpretation  of  certain  important  aspects  of  the 
problem;  these  are  discussed  at  some  length  in  Section  IV..  In  particu¬ 
lar,  the  reason  for  the  introduction  of  logarithmic  terms  into  the  so¬ 
lution,  which  remains  somewhat  of  a  mystery  in  the  artificial-parameter 
analysis,  is  explained  simply  by  the  present  analysis. 

There  is  one  result  of  the  Childress  analysis  which  is  of  special 
significance  in  the  formulation  of  the  far-wake  problem  within  the  con¬ 
text  of  boundary- layer  theory.  The  pressure  gradient  which  appears  in 
the  boundary- layer  equations  would  have  to  come  from  the  solution  out¬ 
side  the  wake.  For  the  "needle,"  as  for  the  corresponding  flat -plate 
problem,  the  pressure  gradient  may  be  taken  as  zero.  Childress shows 

that  the  lowest-order  contribution  to  the  pressure  gradient  begins  af- 
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fecting  the  velocity  distribution  u  at  the  x  term,  but  without  know¬ 
ing  more  about  the  initial  profile  of  u  than  the  total  drag,  this  term 

-2 

is  indeterminate.  Hence,  up  to  but  exclusive  of  the  term  of  order  x  , 
which  includes  all  the  determinate  terms  in  the  expansion,  the  asymp¬ 
totic  expansion  for  u  presented  in  this  Memorandum  is  valid  not  only 
for  the  "needle,"  but  for  any  finite,  closed,  axisymoetric  body. 
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II.  BASIC  EQUATIONS  AND  ASYMPTOTIC  SOLUTION 


We  introduce  a  cylindrical  coordinate  system  (x^,  r^),  where  x^ 
is  the  axial  distance  measured  from  the  base  of  the  cylinder  and 
is  the  radial  distance  measured  from  the  axis  of  symmetry.  The  corre¬ 
sponding  velocities  will  be  denoted  by  u^  and  v^,  respectively. 

In  this  coordinate  system,  the  continuity  and  momentum  equations 
without  pressure  gradient  are 


3^  (rlul>  +  3^  <Vl> 


du^  du^ 

U1  dx^  +  V1  dr^ 


rl  Brl  \  1  &rl 


These  are  to  be  solved  subject  to  the  boundary  conditions 


ux  -  U 


-  0 


and  in  the  limit  of  x.  -* 

(7)  1 

Childress,  in  his  asymptotic  analysis  of  the  Navier-Stokes 
equations  for  the  flow  at  large  distances  from  a  finite  axi  symmetric 
body,  has  shown  that  the  pressure  in  the  inner  viscous  wake  has  the 
following  asymptotic  form: 


a  1  2  .  3  .  x 

4^  72  e  +  °Ce  Xn  e) 
x 


where  x  *  ex/L,  P  is  the  uniform  free-streara  pressure,  L  is  a  charac¬ 
teristic  length,  a  is  a  constant  related  to  the  dimensionless  drag, 
and  € ,  an  "artificial  parameter,"  is  the  ratio  of  the  characteristic 
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length  to  the  length  of  an  extraneous  standard  of  measurement.  Physi¬ 
cally,  the  leading  term  in  this  expansion  is  the  pressure  induced  by 
a  potential  source,  representing  the  effect  of  the  displacement  thick¬ 
ness  on  the  outer  flow.  We  see  that  far  from  any  finite  axisymmetric 

-2 

body  the  pressure,  up  to  but  not  including  the  term  of  order  x  ,  is 
constant;  consequently,  to  that  order  the  formulation  of  the  far -wake 
problem  as  given  above  for  zero  pressure  gradient  applies  not  only  to 
the  "needle"  but  to  finite  bodies  as  well. 

Equation  (1)  is  satisfied  by  introducing  a  stream  function  ^  de¬ 
fined  by 


rlUl 


st; 


r  v 
ii 


dx. 


(5a) 


(5b) 


It  is  now  convenient  to  introduce  a  dimensionless  stream  function 
f(C,xj)  such  that 

♦  l  -  vx1f(C,x1)  (6) 

where 


e 


(7) 


Then  the  nor.dimensional  velocity  components  u  and  v  can  be  written: 


u 


1 


U 


u  = 


(8) 
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v 


U 


(£  -  Cfc  +  .jf^) 


(9) 


(The  subscripts  are  derivatives.) 

In  terms  of  these  new  independent  and  dependent  variables,  Eq.  (2) 
becomes 


CfCC c  +  fQC  +  2  ff££  +  2  Xlfx1fCC  "  2  Xlfx1CfC  "  0  (10) 


Now  x^  is  nondimensionalized  by  means  of  the  momentum  thickness  0^,  de¬ 
fined  by 


tt0^U^  =  2tt  J  (U  -  u^)u^r^  dr^ 

o 


(ID 


or 


4ux  °° 

U-u)«aC 


(12) 


Here  0^  is  a  constant,  related  to  the  drag  of  the  body,  D^,  by 


2,2 


Dj_  -  Tip^  0X 


(13) 


We  now  set 


(14) 


The  differential  equation  Eq.  (10),  however,  is  unaffected  by  this 
change,  since  it  is  homogeneous  in  x^. 
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For  the  asymptotic  solution,  we  assume  that  f(£,x)  has  an  ex¬ 
pansion  for  large  values  of  x  of  the  form 


£  f(C) 

f(c,*)  -  l  -V 


(15) 


n=0 


Substituting  into  Eq.  (10)  yields 


00  OO 


®  00  .. 
nf  _  f" 


oo  *  oo  / 

nf  r-,  f 


n  x  _  x 

n»0  n*0 


„  x  _  x 
n=0  n=0 


rr=0  x  n*0X 


n»0  x  n-0  X 


(16) 


Equating  the  coefficients  of  powers  of  x11  to  zero  yields  the  following 
set  of  ordinary  differential  equations: 


£f '"+  i"  +  ^  f  i" 
o  o  2  o  o 


(17) 


Cfl'+f'l+2fof'l  +  2  £o£l  -  0 


(18) 


Cf2"+  t\  +  \  it\  -  \  f'0t2  +  £'£'  +  \  t[2  -  0 


(19) 


££3'+  f'3+J  £0£'  +  f  £'£' 


-  f  f' 

2  r2rl 


f#f  +~f/f/ 
or3  2  rlr2 


0  (20) 


Since 


.  r  fn(c) 

♦  i  -  vx  )  - - 

io  <vv 


we  note  that  the  boundary  condition*,  Eq.  (3),  require 
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£n(0)  -  0 


-  0  as  C  -  0 


Boundary-condition  Eq.  (4)  Implies  that 


00 


I 

n»0 


£'(« 


—  2  as 


The  condition 


-*  U  as  -*  « 


independent  of  £ 


must  also  be  satisfied,  and  this  indicates  that 


fQ(£)  -  2  as  £  -»  ® 


Equation  (23)  then  reduces  to  the  condition 


f  '(G)  -  0  as  £  - 


n  £  1 


SOLUTION  FOR  f  (O 


XT 


The  differential  equation  is  Eq.  (17): 


III  ,,  1  // 

£f  +  f"  +  ±  f  f  =  0 
o  o  2  o  o 


with  boundary  conditions 


--  ^rr 


(21) 

(22) 

(23) 

(24) 

(25) 

(25a) 

(26) 
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fo(0)  -  0 

V?f*(0  -*  0  as  C  "*  0  ► 


f'(0  -  2  as 


C  -  *  J 


(27) 


A  solution  satisfying  the  differential  equation  and  all  the  boundary 
conditions  is  obvious  by  inspection  and  is 


f0(C)  -  2C 


(28) 


SOLUTION  FOR  f  (Q 


Substituting  Eq.  (28)  into  Eq.  (18)  for  f^(C),  we  obtain 


CfJ>  <C  +  1>£1  +  fi 


(29) 


the  solution  of  which  must  satisfy  the  boundary  conditions 


fL(0)  *  0 


y/Qf[(0  -*  0  as  C  -*  0 

f^(C)  0  as  C  -*  00 


(30) 


If  we  set 


f{(0  -  e'Cgl(C) 


then  g^(Q  satisfies  the  equation 
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Q&l  -  (C  -  1)81  -  0 


whose  solution  is 


*!<»  -  cJ^dC  +  C! 


where  c  and  are  arbitrary  constants. 
Therefore, 


f'(C)  -  ce"C  dC  +  cte'C 

It  is  easily  seen  that  the  second  boundary  condition,  Eq.  (30),  cannot 
be  satisfied  unless  c  =  0;  hence 

f'(Q  -  Cle“c  (31) 

This  automatically  satisfies  the  third  boundary  condition.  Integra¬ 
ting  and  applying  the  remaining  boundary  condition  yields 


f^C)  -  CjU  -  e‘C) 


(32) 


The  constant  c^  is  directly  related  to  the  drag  of  the  body 
through  9^.  In  particular,  since  9^  is  independent  of  x,  we  can  write 
Eq.  (12)  in  the  limit  x  -*  00 , 


9 


2 

1 


lira  f"  .. 

x.-»  U  J  2x  dt 

l  o 


-4u91 

U 


dC 


and  obtain 
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ue. 

2v 


(33) 


SOLUTION  FOR  f  (O 

Using  Eqs.  (28)  and  (31),  we  can  write  Eq.  (19)  as 


Cfj  +  (C  +  l)f"  +  2f2'  +  \  c2  e'2C  «  0 


(34) 


subject  to  the  boundary  conditions 


f2(0)  -  0 


vffJ(C)  -  0  as  £  -  0 


fo(C)  -  o  as  C  “* 


00 


(35) 


Let 


f2(C) 


(C  -  l)e’Cg2(C) 


(36) 


Then  g2  satisfies 


C(C  -  Dg2  +  C2C  -  (C  -  1)  ]g2 


1  2  -C 

2  C1  e 


This  is  easily  integrated  and  yields 


g2(C)  -  4 — - - o  C(C 

4  c(c  -  ir 


•j)e  +  c'] 


If  we  integrate  this  result  and  substitute  in  Eq.  (36),  we  obtain 
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f2'(0 


-t 

e 


t 


+  c 


2 


where  A  and  are  arbitrary 
As  £  “*  0 


constants. 


rC  e-t 

— —  dt  ~  in  £ 
JA  C 


rC  ec 

— r  dt  ~  Xn  £ 

JA  C 

which  means  that  f'(£)  “  as  4n  £  when  £  -*  0-  but  this  would  indi 

cate  that  u  is  infinite  along  the  axis,  ^  -  0.  To  avoid  this,  we 
must  choose 


c 


2 


(37) 


The  boundary  condition  f^CC)  -*  0  as  £  -*  00  is  satisfied,  since 
-C  eC 

£e  ■=-  dt  -*  1  as  £  -*  ® 

JA  C 

One  can  also  verify  that  f"  constant  as  £  -*  0,  so  that 

Vcf2(0  -  0  as  £  -  0 


w  * 


s.  - — — y 
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as  required  by  the  middle  boundary  condition  of  Eq.  (35).  This  com¬ 
pletes  the  determination  of  f^CO*  which  may  be  written  as 


f2'(0 


(C  -  l)e‘c 

\c  t  -t 

1  - — dt  +  2c  -  2 

-  e‘2C  -  l\ 

l 

V  " 

0 

> 

(38) 


(The  boundary  condition  *  0  is  obviously  superfluous  as  far  as 

the  determination  of  f'  is  concerned.)  Note  that  after  satisfying  all 
the  boundary  conditions,  an  arbitrary  constant  c  still  remains  in  the 
equation  for  f^Q* 

To  second  order  then,  the  velocity  on  the  x^  axis  is  given  by 


uo(x) 


0  0  0 


where  c  is  an  undetermined  constant. 
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iii.  modification  of  expansion 


The  occurrence  of  an  arbitrary  constant  in  the  expression  for 
f'(Q,  and  consequently  for  uQ(x)  ,  is  not  in  itself  surprising.  In 
fact,  since  the  initial  conditions,  ignored  in  calculating  the  asymp¬ 
totic  expansion,  influence  the  far-wake  solution,  we  must  anticipate 
that  at  some  stage  of  the  expansion  indeterminacy  will  occur Of 
greater  concern  is  the  fact  that  the  second  approximation,  f'(0,  does 
not  approach  zero  exponentially  as  C  "*  “»  as  we  would  expect,  since 


lim 

C-® 


lim 

“  C-* 

A 

(C  -  De 


-c 


We  note  that  f '(C)  goes  to  zero  as  C’2  a»  C  -  *•  A  si">ilar  Phenomenon 
was  found  by  Goldstein^  in  the  two-dimensional  oase^in  the  third 
approximation  (his  solution  proceeded  in  terms  of  x  ”  ,  n  taking  on 

all  integral  values).  Even  if  this  approximation  were  not  rejected 
on  physical  grounds,  it  is  actually  impossible,  as  shown  by  Goldstein, 
to  satisfy  all  the  boundary  conditions  in  the  next  approximation,  so 
the  solution  is  not  even  possible  mathematically.  It  is  very  likely 
we  would  encounter  the  same  difficulty  if  we  attempted  to  solve  for 

the  next  term  in  our  expansion,  f^CO* 

Stewartson(5)  resolved  the  difficulty  for  the  two-dimensional 
flat-plate  case;  he  showed  that  the  asymptotic  expansion  contained  an 
extra  term  0(x"3/2  4n  x)  in  addition  to  the  x  n  terms  (n  =  1,  2,  3) 
contained  in  Goldstein's  original  solution  and  also  that  a  numerical 
factor  occurring  in  the  term  0(x  3^2)  was  indeterminate.  Stewartson 
treated  a  number  of  boundary- layer  problems  involving  asymptotic  ex¬ 
pansions  for  large  x  and  demonstrated  that  these  phenomena  (the  occur 
rence  of  logarithmic  terms  and  indeterminacy  in  the  following  alge¬ 
braic  terms)  are  general  features  of  all  of  them. 
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Brief  ly,  Stewartson's  demonstration  proceeds  as  follows:  In  each 
of  the  problems  considered  at  the  nth  stage  of  the  expansion  procedure , 
the  required  function  satisfies  a  third-order  ordinary  linear  differ- 
ential  equation  with,  say,  7]  as  an  independent  variable  and  n  as  a 
parameter;  two  boundary  conditions  are  given  at  T)  =  0  and  one  at  T|  « 

Ft  each  of  an  infinite  countable  set  of  n's,  not  necessarily  integers, 
there  is  a  complementary  solution  satisfying  all  the  boundary  condi¬ 
tions  at  7]  «=  0  and  at  T]  =  ».  (These  are  called  "eigensolutions."^6,8^) 
It  follows  that  the  factor  multiplying  this  function  cannot  be  deter¬ 
mined;  and  Stewartson  concludes  that  it  probably  depends  on  the  neg¬ 
lected  upstream  boundary  condition,  the  initial  velocity  profile.  When 
a  member  of  this  set  is  an  integer,  then  no  particular  integral  of  the 
differential  equation  is  exponentially  small  for  large  7]  unless  an  ad¬ 
ditional  term,  consisting  of  the  complementary  solution  multiplied  by 
ln  X  °r  ln  in  x»  depending  on  the  problem,  is  added.  The  factor  in 
this  term  is  determined  by  requiring  that  the  particular  integral  be 
exponentially  small  for  large  T|;  the  particular  integral  so  defined 
is,  however,  not  unique. 

In  the  following,  we  shall  demonstrate  the  particular  form  which 
the  ideas  expressed  above  take  in  the  present  problem.  Using  the  same 
general  approach  as  Stewartson  employed  for  the  flat-plate  problem, 
we  shall  indicate  why  the  expansion  of  the  previous  section  broke  down 
and  proceed  to  determine  the  correct  expansion.  In  this  process,  we 
shall  explicitly  determine  the  set  of  n  tor  which  eigenfunctions  exist 
and  the  nature  of  the  eigenfunctions  themselves.  We  shall  see  that 

there  are  significant  differences  between  the  two-dimensional  and  axi- 
symmetric  cases. 

Consider  Eqs .  (1)  and  (2)  again  with  the  boundary  conditions 


uL  -  U 


as  r  • 


xl>  0 


r  -  h 

1  dr. 


r  <39> 


-  0 


at 


-  0 


x1  >  0 


j 


-16- 


and  the  Initial  condition 

«  UfCr^  at  =  0  >  0  (40) 

This  initial  condition  is  introduced  here  only  for  convenience  in 
carrying  out  the  later  discussion,  which  will  lead  to  the  correct  form 
of  expansion.  This  initial  condition  should  not  appear  in  the  actual 
formulation  of  the  asymptotic  problem.  (It  does,  however,  play  a  role 
in  determining  certain  constants  in  the  asymptotic  solution.)  In  gen¬ 
eral,  an  axisymmetric  body  of  nonzero  thickness  will  have  a  separated 
region  immediately  behind  the  base,  and  the  flow  field  is  too  complex 
to  allow  the  computation  of  f(r^).  However,  for  the  purposes  of  our 
later  discussion  it  is  not  necessary  that  the  initial  profile  be  given 
at  the  base  of  the  body;  it  is  only  necessary  that  some  initial  pro¬ 
file  be  specified,  and  we  arbitrarily  place  the  origin  of  the  x^  axis 
at  the  place  where  the  profile  is  given.  Thus,  for  example,  we  can 
think  of  f(r^)  as  being  given  at  some  location  downstream  of  the  sep¬ 
arated  region  if  one  exists.  In  any  case,  whether  or  not  f(r^  can 
readily  be  determined  is  immaterial  to  our  later  arguments.  However, 
we  do  assume  that  f(r^  is  some  sort  of  boundary- layer  solution.  As 
a  corollary  of  this,  it  follows  that  f(r^  -*  1  exponentially  as  r^  ap¬ 
proaches  infinity. 

We  nondimensionalize  by  setting 


(41) 


Equations  (1)  and  (2)  now  become 


h  <ru)  +  I;  (rv)  ' 


u 


i  a  /  du\ 
r  dr  \r  dr  ) 


(42) 


du 

dr 


(43) 
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A  first  approximation  to  u  when  x  is  very 
from  the  Oseen  approximation  and  the  fact  that 
is  constant  (since  the  pressure  is  taken  to  be 
is 


large  can  be  obtained 
the  momentum  thickness 
constant).  The  result 


u 


A  -r  /4x 

-  e 

x 


u. 


(44) 


say,  where  A  is  Riven  in  terms  of  the  nondimens ional  drag  D  by 


D 


4npA 


(45) 


or 


A 


f(l  -  f)r  dr 


(46) 


To  improve  this  approximation  we  substitute  Eq.  (44)  into  Eqs. 
(42)  and  (43),  thus  starting  an  iteration.  To  demonstrate  the  nature 
of  this  process  at  any  stage,  we  set 


u  ■ 


v  ■ 


(47) 


Substituting  into  Eq.  (43),  we  rewrite  it  as 


du  1  d  /  du  \  —  du  —  du 

dx  "  r  dr  dr  j  "  ’  U  dx  "  V  dr 

th 

The  solution  in  the  n  iteration  is  then  obtained  from  the  equation 


— (n-1)  du^n  — (n-1)  du^n  ^ 

u  dx  "  V  dr 


(48) 


’flu 
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(This  iterative  method  of  solution  can  be  shown  to  be  equivalent  to 
the  asymptotic  series-expansion  approach  discussed  in  Section  II.) 
For  n  ■  1,  taking  u^°^  -  0,  Eq.  (48)  reduces  to 


auW 

dx 


0 


(49) 


the  solution  of  which  is 


-(1)  A  -r  /4x 

u'  - e 

x 


(50) 


(already  given  in  Eq.  (44)). 

When  n  =  2,  Eq.  (48)  reads 


dx 


i  a_ 

r  dr 


-(1)  — ^ 

-  “CD  (51) 


dx 


dr 


By  substituting  Eq.  (50)  into  Eq.  (42),  we  find  that 


-r2/4x 

e 


(52) 


Substitution  of  Eqs.  (50)  and  (52)  into  Eq.  (51)  now  yields  as  the 
governing  equation  for  the  second  approximation 


ST<2>  nf  c*<2>\  A2  -r2/2x 
Si  ’  i  Jr  V  9r  '  3  e 

X 


(53) 


Rather  than  determine  the  solution  to  Eq.  (53)  at  this  time,  we 
want  to  consider  the  general  problem  at  any  stage  of  the  iteration. 
In  each  iteration,  we  are  treating  the  equation 

du  _  1  d_  /  du\ 
dx  r  dr  \  dr  / 


F(x,r) 


(54) 
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where 


F(x,r)  -  (1  -  «)  f*  -  v  £  (55) 

the  value  of  F  being  calculated  using  the  u  and  v  determined  from  the 
solution  of  Eq.  (54)  in  the  previous  iteration. 

Equation  (54)  is  of  the  form  of  an  inhomogeneous  heat  equation, 
and  the  determination  of  its  solution  is  a  standard  problem  in  heat 
conduction.  Assuming  f(r)  and  F(x,r)  to  be  even  functions  of  r,  the 
solution  of  Eq.  (54)  satisfying  the  boundary  conditions,  Eqs.  (39)  and 
(40),  is  given  by 


u 


2x 


e 


(r2+r/2) 

4x 


dr' 


+ 


dr 


2(x  -  x') 


e 


(r2+r'2) 

40M7) 


(56) 


where  Iq  is  the  modified  Bessel  function  of  the  first  kind  and  zeroth 
order. 

For  simplicity  and  without  sacrificing  any  of  the  significant 
details,  we  shall  limit  ourselves  to  determining  u  on  the  axis,  r  =  0. 
Setting  r  =  0  in  Eq.  (56)  yields 


u 


r=0 


=  1 


1  -  f  Cr  ') 

2x 


r 

4x 


*  ,  * 
r  dr 


r-'/ 


dr'  rlFCx.V')  4(x-x^) 
Qr  2 (x  -  x')  6 


(57) 
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With  the  general  solution  for  any  iteration  now  given  in  terms 

_(2) 

of  quadratures,  we  return  to  the  determination  of  uv  .  From  Eq.  (53) 
we  note  that  F(x,r)  can  be  written 


where 


F(x,r)  =  ~  e“r  ^2x  =  g(x)F2(T)) 

x 


(58) 


g(x)  =  ~ 


x 


F2(T1) 


-T 


(59) 


Introducing  Eq.  (58)  and  the  new  variable  T|  into  the  second  integral 
of  Eq.  (57),  we  obtain 


u 


1 


00  ^ 

r  1  -  f (r  ) 

J  2x 


e 


dr 


* 


+ 


^  g  (  x  *)  x  *  dx  * 
J  (x  -  x') 


H^COe 

'  o 


2<x  x'>  dll' 


(60) 


Now  we  write 


l'(x) 


xg(x) 


(61) 
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The  contribution  from  F  is  then 


xV 

2  (x-x 


2 

1 


<m' 


Integrating  by  parts,  this  becomes 


(62) 


x 


VF2<TT)e 


xV2 

2(x-x') 


(63) 


(We  cannot  conclude  that  L(0)  is  infinite  from  Eq.  (61),  since  the  it¬ 
erative  solution  based  on  Eq.  (54)  is  valid  only  for  large  x.)  For 
large  x  the  leading  terms  in  Eq.  (63)  are 

-  ^  J  Vf2OT)  dll'  +  ^  2n  X  [  (l  -  2^)T|'F2(T|')  dll'  (64) 

O  X  o  \  / 


Following  the  method  of  Stewartson,  two  questions  can  now  be 
raised.  First,  does  the  initial  profile  affect  u  only  through  the 
constant  A,  given  in  terms  of  f(r)  by  Eq.  (46)?  Second,  is  it  pos¬ 
sible  to  find  a  second-order  approximation  which  has  the  correct  asymp¬ 
totic  behavior  as  r  -*  ®,  thus  eliminating  the  difficulty  encountered 
in  Section  II;  and  if  so,  what  is  the  correct  second  approximation? 

The  first  question  may  now  be  answered  on  the  basis  of  expression 
(64).  The  first  term  represents  a  modification  to  the  term  0(l/x)  of 
the  expansion  for  large  x;  this  occurs  because  the  contribution  from 
f(r)  in  Eq.  (60)  to  the  0(l/x)  term  is 


Sj  a-f<r'»*'dr' 


(65) 
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which  differs  from  Eq.  (46).  There  will  be  a  contribution  to  expres¬ 
sion  (65)  at  each  stage  until  the  coefficient  of  1/x  is  Eq.  (4  6) .  In 
Eq.  (57)  the  general  term  of  the  contribution  from  f(r)  to  u(x,0)  is 


("1)  P  /.  f/  *\\  '2rtfl  ,  *  trr\ 

— 4 — *■—  (1  -  f (r  ))  r  dr  (66) 

2x(4x)  n.'  J  o 

and  in  a  similar  way  this  will  be  modified  at  each  stage  of  the  itera¬ 
tion.  Thus,  we  can  expect  terms  in  the  expansion  of  order  x  n  ^  to 
have  an  indeterminate  numerical  factor  depending  on  f(r).  The  term 
with  n  =  0  is  completely  determined  because  of  the  momentum  integral, 
but  this  is  the  only  such  integral  known. 

Expression  (64)  enables  us  to  answer  the  second  question  also. 

The  second  term  of  expression  (64)  shows  that  the  asymptotic  expansion 

2 

contains  a  terra  of  the  form  in  x/x  ;  it  was  the  absence  of  this  term 
in  the  expansion  in  x  n,  considered  in  Section  II,  which  caused  the 
difficulties  encountered  in  the  second-order  approximation. 

Let  us  now  return  to  Eq.  (53)  and  assume  that 

-(2)  _  injc  G(T|)  +  i_  H(1))  (67) 

X  X 


Substitution  in  Eq.  (53)  yields  the 
TIG"  +  ( T) 2  +  1)G  1  +  4T)G 

T)H"  +  (T\2  +  1)H 7  +  4T|H 


equations 


=  0 


-  2T|G  -  2A 


(68) 

(69) 


The  boundary  conditions  are 
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G7(0)  -  H7(0)  -  0 


and 


(70) 


G,H  -*  0  exponentially  as  T)  -*  ®  J 


The  solution  of  Eq.  (68)  satisfying  the  boundary  conditions,  Eq. 
(70),  is 


G(T» 


(71) 


Substituting  this  result  into  Eq.  (69)  and  changing  to  the  new  inde¬ 
pendent  variable  £  defined  by 


c 


(72) 


(this  £  is  the  nondimens ional  fora  of  the  £  used  in  Section  II),  we  ob¬ 
tain 


£H"  +  (£  +  1)H 7  +  2H  -  C(£  -  l)e'C  -  A2  e‘2C  (73) 

where  primes  now  denote  differentiation  with  respect  to  £.  This  can 
be  solved  by  setting 


H(£) 


(£  -  l)e_CK(£) 


Here,  K(£)  satisfies  the  equation 


K"  +  ^  -  l)2]  K,  _  C(£  - 


1)  -  A2e~£ 


£(C  -  1) 


za  - 1) 


(74) 
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A  first  integral  of  this  equation  is 


K'  e"C  C(C  "  l)2  =  JC  e"C  (C 

"o 


(C  “  1) 


[c(C  -  1)  -  A2  e"CJ  dC 


If  H(T1)  is  to  tend  exponentially  to  zero  as  7]  -*«(£-♦  ®)  ,  then  we 
must  require 

J  e'^  (X  -  1)  [c(C  -  1)  -  A2e'^]  d£  «  0  G 


This  condition  determines  C  to  be 


1A2 

4  A 


Integrating  Eq.  (75)  once  more  and  substituting  into  Eq.  (74), 
we  obtain  the  H(Q  satisfying  the  differential  equation,  Eq.  (69), 
and  boundary  conditions,  Eq.  (70): 


H(0 


dt  +  e 


+  4  [« - 


l)e"£  In  l  -  2e" 


where  a  is  an  undetermined  constant. 

Thus,  G  is  determined  uniquely  but  H  is  not;  this  is  not  sur¬ 
prising,  since,  as  we  saw  earlier,  H  probably  depends  on  f(r). 

The  solution  for  u  to  the  present  order  of  approximation  is 
given  by 


«  -  1  -  A.'r  /4x  +  cCo  +  1_  H(o 
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2  o 

where  £  =  T|  /2  =  r  /4x  and  G(£)  and  H(£)  are  given  by  Eqs.  (71)  and 
(78),  respectively.  Thus,  we  can  also  write  the  solution  as 

u  =  1  -  A  e-^  -  I  A2(C  -  l)e’^  I?-*  +  ~  H(£)  (80) 

x  x 

On  the  axis  C  =  r  **  0,  this  reduces  to 

u(x,o)  =  1  -  %  a2  +  L  (81) 

X  x 


where  8  is  an  undetermined  constant. 

Let  us  now  consider  what  happens  in  the  next  iteration,  the  de- 
termination  of  uv  '  .  Equation  (48)  now  reads 


1  3 


dx  r  dr 


(■#) 


-(2)  du(2)  -(2)  du(2) 

U  dx  "  V  dr 


(82) 


For  convenience  in  calculating  the  right-hand  side,  we  introduce  1 
r/s/Sx  as  an  independent  variable;  this  allows  us  to  write 


-  u 


77(2)  du(2)  -(2)  du(2) 


dx 


dr 


"(2)  dui^  (ll  -(2)  _ 

**  +\2xU 


du<2> 

dTl 

(83) 


In  terms  of  T|,  Eq.  (80)  becomes 


u  =  !  -  f  *L5  +  i-  „„)  (84) 

'  /  XX 


It  will  be  convenient  for  our  purposes  if  we  represent  Eq.  (84)  by 


»,<H>  In  x  M1)) 

u  -  1  -  -V  -  V>  *¥  -  -V 

x  x 


1  +  u 


(85) 
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Integrating  the  continuity  equation,  Eq.  (42),  and  using  Eq.  (85), 
we  find 


v 


-2 

nyjx 


1  -2  .  r  in  x  .  G3^ 

u  V  +  g2(11)  ~t  +  — r 

X  X 


(86) 


where  G0  and  G„  are  given  in  terms  of  D.  and  D  but  their  explicit 

L  j  L  j 

representation  is  not  necessary  at  this  time. 

If  we  now  substitute  Eqs.  (85)  and  (86)  into  the  right-hand  side 
of  Eq.  (83),  this  equation  becomes 


-(2)  du(2)  -(2)  8u(2)  1  „2  .  in  x 

-  u  “ST  -  v  -5T  ‘  “  D1  +  — 


3D1D2  •  Dl'l  T  -  2  T'D2j 


X 


3did3  +  did2  -  Di  (-  5 ltD: 


3  Ti 


+  . . 


(87) 


The  right-hand  side  of  Eq.  (87)  is  now  substituted  for  the  right-hand 

side  of  Eq.  (82)  and  is  therefore  the  inhomogeneous  term  for  the  third 
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iteration.  The  term  D^/x  is  exactly  the  one  treated  earlier  in  the 
second  iteration  (see  Eq.  (58)),  and  since  F(x,r)  enters  Eq.  (57)  (or 
Eq.  (60))  linearly,  it  will  contribute  only  the  terms  obtained  in  the 
second  iteration.  Hence,  we  shall  limit  our  attention  to  the  remain¬ 
ing  two  terms  on  the  right-hand  side  of  Eq.  (87). 

We  now  let 


F(x,r)  -  g1(x)F;J(Tl)  +  g2(x)F4(T]) 


(88) 


gL(x) 


In  x 

4 

x 


g2Cx> 


where 
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and  and  are  the  coefficients  in  Eq.  (87),  multiplying  and  g2> 
respectively.  Correspondingly,  we  define 


so  for  large  x 


L^(x)  «  xg^x)  = 

x 


L2;(x)  -  xg2(x)  =  —■ 

X 


„  ,  x  1  In  x  1 

ll(x)  ~  -  2  — 2 - 2 

x  4x 


L9(x) - r 

1  2yT 


Now  consider  the  second  term  of  Expression  (63) ,  giving  the  con¬ 
tribution  from  F(x,r)  to  u;  in  particular,  we  take  the  sum  of  two  such 
terms,  L^(x)  and  appearing  in  one  and  L2(x)  and  F^CH)  in  the 


other: 


jc  L  ,(x')  dx  '  [  _^2 

-  f  7 - 75  W  1  ' 

o  (x  -  x  )  o  L  J 


XV  n'F  r-n')e~  2(x"X  } 
2 (x  -  x  1  V  1  )e 


nX  L  (x')  dx'  ® 

- 1  -f - 72  !  *n'  ] 

o  (x  -  X  )  o  L 


1  -  1  *4<0.  ^  > 


For  large  x  the  leading  terms  are 


h\\  Li(x'>  dx\l 

X  o  o 


r  t  -  4) 


Tl'F3(7)')  dfl' 
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+  I 

x 


f  Li<x'>*' «*'  J  (2  - 

0  o  \ 


27] 


'2  +  al 

4 


^F3(71")  dTj" 


+  J*  l2(x")  dx'  J  (l  -  ^tTf^tT)  dll' 


+  1 
x 


j  L2(x')x"  dx' 


^2  -  27i '2  + 


(90) 


According  to  Eq.  (89) , 


r 


T  ,  '\  ,  *  in  x  ,  3 
Li(x  )  dx  +  ^ 


r 


L2(x')  dx'~ij 


r 


L..(x')x'  dx' - 7  (2n  x)^  -  7  -Cn  x 

1  44 


L2(x')x'  dx' 


j  in 


Substitution  of  these  integrals  into  expression  (90)  allows  us  to  write 
it  as 
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2  -  27)  + 


'2  4.1^ 


1  Tl  'f,  (T1 ')  dn'  + 


in  x 
2x3 


-  ^tTf^tThtT 


1 

2 


"  (2  -  2^1  '2  <■  VVF3(T1")dTl'  -  J  (2  -  2T1/2  +  ^)t]/F4(11')  d1! ' 


2x 


I  J^1  -  )vF3oi')dV  +  J^I  -  d-n' 


(91) 


(3) 


Thus  we  see  that  in  the  solution 'for  u  ,  terms  of  order 

2  3  3  3 

(in  x)  /x  ,  in  x/x  ,  and  l/x  appear. 

To  gain  some  insight  into  why  terms  involving  powers  of  the  natu¬ 
ral  logarithm  should  appear  in  the  expansion,  let  us  return  to  Eq.  (48). 
The  complementary  solution  at  any  stage  of  the  iteration  is  determined 
as  the  solution  of  Eq.  (48)  with  the  right-hand  side  set  equal  to  zero. 
In  terms  of  the  variables  x  and  T]  =  rA/2x,  the  complementary  function 
satisfies  the  equation 


1  a  L  su(n) 

dx  2  x  ST]  j  2xTl  37]  \ 11  3T) 


(92) 


If  we  now  substitute 


-(n)  V™ 

u  ■  - 


X 


n 


(93) 


we  obtain  the  following  equation  for  Q  (T|)  : 

n 


T1Q"  +  01  +  1)Q  '  +  2nT)Q  -  0 
n  n  n 


(94) 


Now  let 


Qn(T1)  "  Mn(T1)e 


-Tl2/2 


(95) 
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This  yields 

Mn  +  (^  •  ^  Mn  +  2(n  ‘  1)Mn  '  0 

2 

Finally,  we  again  introduce  the  independent  variable  £  =  T)  /2  in  terms 
of  which  the  equation  for  becomes 

GM"(C)  +  (1  -  OM'(C)  +  (n  -  1)M  (0  -  0  (96) 

n  n  n 

The  solution  of  this  is 

V«  -  L„-i(C) 


where  L  (£)  is  the  Laguerre  polynomial  of  order  n.  Hence 
n 


Q  (Tl)  =  L  (\)  e_T1  /2 

n  n-l\  2  I 


The  boundary  conditions  to  be  satisfied  by  Qn01)  are 

0^(0)  -  0  Q  (T])  -♦  0  exponentially  as  T]  -»  oo 

n  n 


The  second  boundary  condition  is  immediately  seen  to  be  satisfied. 
For  the  first  condition,  we  note  that 


q'CH) 


e'^2  -  TIL 


-T|2/2 


so  Q^(0)  ■  0,  and  this  condition  is  also  satisfied  for  any  n.  Hence, 
the  complementary  functions  for  all  n  exactly  satisfy  the  boundary 
conditions  at  T)  ■  0  and  7]  «=  «. 

The  corresponding  inhomogeneous  equation  looks  like 
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T1R"  +  (T|2  +  1)R^  +  2nT]Rn  =  F(7))  (98) 

where  F(T))  is  a  known  function  of  T|,  perhaps  involving  arbitrary  con¬ 
stants.  If  we  let 


Rn(T»  “  S^TDe'11  12  (99) 

2 

and  also  introduce  the  variable  Q  =  7]  /2  again,  the  equation  becomes 

cs;«)  +  (1  -  0S'<0  +  (n  -  l)Sn(0  -  (100) 

In  particular,  we  consider  the  third  approximation,  n  *  3.  The 
equation  is  then 


CS"  +  (1  -  QS'  +  2S3 


v.LQ£ 

iM 


(101) 


If  we  set 


s3(0  =  l2(Qt(0  -  (C2  -  4C  +  2)T(0  (102) 


(L2(Q  is  the  Laguerre  polynomial  of  ore  *r  2),  then  we  obtain 


T 


// 


f-C3  ±  9£2  -  14 C  +  2^, 
\  C( c2  -  4£  +  2)  / 


F  (OeC 

2qVTq  (C2  -  4C  +  2) 


The  first  integral  of  this  equation  is  given  by 


e"CC(C2  -  4C  +  2)2  t' 


^  F(Q(C2  -  4C  +  2) 

2/2Q 


dC 


(103) 


Thus,  if  R3  is  to  tend  to  zero  exponentially  as  Q  or  Tj  -*  we  must 
require 
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r  HQU2  -J&+21  iC  .  o  (104) 

J0  2  J2C 


3  —(3) 

If  R  (TO  represents  the  coefficient  of  the  term  Xn  x/x  in  uv 


that  is, 


u(3)  *  ^  R  (TO  +  term  of  0(x'3) 

x 


then  the  F(Q  to  be  entered  in  Eq.  (98)  and  the  following  equations 

4 

is  the  coefficient  of  Xn  x/x  in  Eq.  (87).  The  explicit  representa¬ 
tion  of  this  in  terms  of  T)  is 


3D1D2 


Hence 


F(C)  -  A3(C  -  l)e‘2C  (105) 

If  we  substitute  this  expression  into  the  integral  in  Eq.  (104),  the 
integral  becomes 


AL  f  (C  -  D(C2  -  4 C  +  2)e~2^ 

2j2  °  A 


3 

The  value  of  this  integral  is  (-89/216)A  ,  and  hence  the  condition 
that  R_  -*  0  exponentially  as  11  -•  ®  is  not  satisfied. 

3  2  3 

This  now  explains  the  appearance  of  the  (Xn  x)  /x  term  in  the 

solution  for  u^ .  This  term  must  be  introduced  to  insure  that  the 

3 

particular  integral  corresponding  to  Xn  x/x  be  exponentially  small 

2  3 

when  T)  is  large.  This  (Xn  x)  /x  term  consists  only  of  the  complemen 
tary  function  (solution  of  Eq.  (94)),  its  numerical  factor  being  de¬ 
termined  from  the  condition  on  the  particular  integral  corresponding 


,  _  ■»  » 


r—  r  yJFT 
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3 

to  in  x/x  .  The  factor  multiplying  the  complementary  function  corre- 

3 

sponding  to  in  x/x  is  determined  by  the  condition  on  the  particular 

3 

integral  corresponding  to  1/x  .  However,  there  is  no  way  of  obtaining 

the  numerical  factor  in  the  complementary  function  corresponding  to 

3 

1/x  in  this  way,  and  as  we  have  already  seen,  it  probably  depends  on 
the  initial  profile,  f(r). 

Generally,  difficulty  arises  each  time  a  complementary  solution 

to  Eq.  (48)  of  the  form  x  nG  (Tl)  satisfies  the  conditions  G*(0)  ■  0 

n  n 

and  -*  0  exponentially  as  T)  -•  ®.  We  have  already  shown  that  this 

occurs  for  all  integers  n,  and  thus  at  each  step  of  the  solution  it  may 
be  necessary  to  introduce  an  extra  factor  in  x  into  the  series.  At 
the  same  time,  the  x  n  term  will  contain  an  indeterminate  numerical 
factor  depending  in  some  way  on  the  initial  profile  of  u. 
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IV,  DISCUSSION 

A  number  of  differences  may  be  pointed  out  between  the  two-dimen¬ 
sional  and  axisymmetric  results.  Stewartson^  shows  that  the  veloc¬ 
ity  at  y  =  0  for  the  flat-plate  asymptotic  far-wake  solution  is 


u(x  ,0) 


where  or  is  a  known  constant  related  to  the  drag  of  the  plate  and  0  is 

an  unknown  constant.  (Stewartson  gives  a  coefficient  of  the  term 

3 /2 

in  x/x  which  is  twice  as  great  as  the  one  given  above;  this  was  due 
to  a  slight  error  in  one  of  his  equations  that  was  corrected  by  Crane.  ^ 
We  have  found  the  corresponding  axisymmetric  solution  to  be  Eq.  (81): 


u(x ,0)  =  1  -|  +  \  A2  *L2i  +  ^+o 


where  A  is  known  in  terms  of  the  drag  of  the  body  and  B  is  unknown. 

We  note  that  the  indeterminacy  in  the  asymptotic  expansion  appears 

earlier  in  the  axisymmetric  than  in  the  two-dimensional  case,  occurring 

in  the  fourth  term  in  the  former  and  in  the  fifth  in  the  latter.  Also, 

in  the  two-dimensional  case,  eigenfunctions  occur  for  terms  involving 

x  n  only  when  n  is  an  odd  integer,  so  a  factor  £n  x  may  be  introduced 

at  alternate  steps  of  the  series  expansion.  In  the  axisymmetric  case, 

on  the  other  hand,  eigenfunctions  exist  for  each  n,  so  these  additional 

logarithmic  factors  may  occur  at  each  step. 

As  indicated  in  Section  II,  the  pressure  in  the  wake  is  constant 
_2 

to  order  x  ;  consequently,  Eqs.  (80)  and  (81)  up  to  terms  of  order 
-2 

x  in  x  represent  the  asymptotic  far-wake  solution  at  large  Reynolds 

_2 

number  for  any  finite  axisymmetric  body.  The  term  of  order  x  ,  which 
is  the  first  term  in  the  asymptotic  expansion  affected  by  a  nonconstant 
pressure,  is  already  indeterminate  due  to  the  existence  of  an  eigenso- 
lution;  thus,  the  determinate  part  of  the  solution  is  unaffected  by  the 
existence  of  a  nonzero  pressure  gradient. 
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It  is  of  some  interest  to  compare  the  approach  taken  here  with 
that  of  Chang^^  and  Childress^^  in  their  analyses  of  the  far  flow 
field  for  arbitrary  Reynolds  number.  As  indicated  earlier,  these  au¬ 
thors  introduce  an  artificial  parameter  e,  which  is  the  ratio  of  a 
characteristic  length  to  an  artificial  length,  and  recast  the  problem 
of  the  flow  at  large  distances  as  a  perturbation  for  small  values  of 
e;  in  other  words,  the  coordinate-type  expansion  is  replaced  by  a  pa¬ 
rameter-type  expansion.  They  then  construct  two  expansions  for  the 
velocity  and  two  for  the  pressure,  an  outer  expansion  valid  outside 
the  wake  region  and  an  inner  expansion  valid  within  the  wake  region. 

The  domain  of  validity  of  these  two  expansions  overlaps ,  and  hence  an 
expansion  uniformly  valid  at  large  distances  may  be  constructed  by  a 
combination  of  outer  and  inner  expansions. 

Apart  from  the  usual  principles  for  matching  inner  and  outer  ex¬ 
pansions,  Chang  and  Childress  employ  two  additional,  independent  prin¬ 
ciples  for  finding  the  form  of  the  expansions  and  for  eliminating  cer¬ 
tain  apparent  indeterminacies.  These  are  (1)  the  principle  of  elimi- 
nability,  and  (2)  the  principle  of  rapid  (transcendental)  decay  of  vor- 
ticity.  The  first  of  these  requires  that  the  artificial  parameter  e 
be  eliminable  from  the  expansions.  The  second  principle,  that  for  a 
finite  or  semi-infinite  solid  in  a  uniform  stream  the  vorticity  decays 
at  an  exponential  rate  with  distance  outside  the  wake  (or  at  least 
faster  than  any  power  of  the  distance)  ,  is  assumed  to  hold  for  solu¬ 
tions  of  the  Navier-Stokes  equations;  all  the  known  exact  and  linear¬ 
ized  solutions  support  this  assumption  (see  Ref.  6).  The  application 
of  the  eliminability  principle  leads  to  the  concept  of  "switchback 
terms"  which  must  be  introduced  in  order  that  e  be  eliminable  from  the 
resulting  expansions.  In  the  Chang  and  Childress  analyses,  the  effect 
of  switchback  is  to  introduce  logarithmic  terms;  thus,  for  example,  in 
the  axisymmetric  solution  of  Childress,  terms  e^(in  e)^  (i  and  j  are 
integers  satisfying  j  £  i  -  1)  are  introduced  into  expansions  which 
originally  contained  only  integral  powers  of  e.  Childress  does  not 
attempt  to  explain  precisely  the  reason  for  switchback  but  notes  that 
it  is  a  nonlinear  phenomenon. 


One  of  the  essential  differences  between  the  artificial-parameter 
analyses  of  Chang  and  Childress  and  the  present  analysis,  based  on  a 
coordinate-type  perturbation,  is  the  reason  for  the  introduction  of 
the  logarithmic  terms.  In  the  present  analysis  these  terms  are  intro¬ 
duced  in  order  that  the  vorticity  decay  exponentially  at  infinity.  In 
the  artificial-parameter  analysis,  as  indicated  above,  the  switchback 
terms  arise  from  the  application  of  the  principle  of  eliminability  of 
the  artificial  parameter;  the  principle  of  rapid  decay  of  vorticity 
is  put  forth  as  a  completely  independent  principle.  Thus,  the  under¬ 
lying  reason  for  the  occurrence  of  the  switchback  terms,  which  is  not 
fully  explained  in  the  Chang  and  Childress  analyses,  would  seem  to 
stem  from  only  one  basic  principle--rapid  decay  of  vorticity. 

It  might  be  added  that  after  listing  certain  of  the  advantages  of 
the  artificial-parameter  approach,  Chang  indicates  that  a  direct  coordi¬ 
nate-type  procedure  would  seem  preferable  for  a  final  analysis. 
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